Differential expression of ASIP transcripts reveals genetic mechanism underpinning black-tail independence from body plumage in yellow-bodied chickens.

The genetic foundation of chicken body plumage color has been extensively studied. However, little attention has been paid to the inheritance patterns and molecular mechanisms underlying the formation of distal feather colors (tail and wingtip). Differences in these colors are common; for example, the Chinese Huiyang Beard chicken has black tail feathers, but yellow body plumage. Here, the hybrid offspring of Huiyang Beard and White Leghorn chickens were used to study the inheritance patterns of tail-feather color. The expression levels of pigment genes in differently colored feather follicles were analyzed using quantitative real-time PCR. The results showed that genetic regulation of tail-feather color was independent of body-plumage color. The Dominant White locus inhibited eumelanin synthesis in tail feathers without affecting the formation of yellow body plumage, whereas the Silver locus had the opposite effect. The expression of agouti signaling protein (ASIP) gene class 1 transcripts was significantly lower in black tail-feather follicles than in yellow body follicles, whereas tyrosinase-related protein 1 (TYRP1) gene expression was significantly higher in black tail feathers. These differentially expressed genes were confirmed to exert an effect on eumelanin and pheomelanin formation in feathers, thus influencing the regulation of chicken tail-feather color. In conclusion, this study lays the foundation for further research on the genetic mechanisms of regional differences in feather color, contributing to a better understanding of plumage pigmentation in chickens.

Role of melanocortin system in the locomotor activity rhythms and melatonin secretion as revealed by agouti-signalling protein (asip1) overexpression in zebrafish.

Temporal signals such as light and temperature cycles profoundly modulate animal physiology and behaviour. Via endogenous timing mechanisms which are regulated by these signals, organisms can anticipate cyclic environmental changes and thereby enhance their fitness. The pineal gland in fish, through the secretion of melatonin, appears to play a critical role in the circadian system, most likely acting as an element of the circadian clock system. An important output of this circadian clock is the locomotor activity circadian rhythm which is adapted to the photoperiod and thus determines whether animals are diurnal or nocturnal. By using a genetically modified zebrafish strain known as Tg (Xla.Eef1a1:Cau.asip1)iim04, which expresses a higher level of the agouti signalling protein 1 (Asip1), an endogenous antagonist of the melanocortin system, we observed a complete disruption of locomotor activity patterns, which correlates with the ablation of the melatonin daily rhythm. Consistent with this, in vitro experiments also demonstrated that Asip1 inhibits melatonin secretion from the zebrafish pineal gland, most likely through the melanocortin receptors expressed in this gland. Asip1 overexpression also disrupted the expression of core clock genes, including per1a and clock1a, thus blunting circadian oscillation. Collectively, these results implicate the melanocortin system as playing an important role in modulating pineal physiology and, therefore, circadian organisation in zebrafish.

The MFSD12 p.Tyr182His common variant is sufficient to alter mouse agouti coat color.

MFSD12 functions as a transmembrane protein required for import of cysteine into melanosomes and lysosomes. The MFSD12 locus has been associated with phenotypic variation in skin color across African, Latin American, and East Asian populations. The frequency of a particular MFSD12 coding variant, rs2240751 (MAF = 0.08), has been reported to correlate with solar radiation and occur at highest frequency in Peruvian (PEL MAF = 0.48) and Han Chinese (CHB MAF = 0.40) populations, suggesting it could be causative for associated phenotypic variation in skin color. We have generated a mouse knock-in allele, Mfsd12Y182H , to model the human missense p.Tyr182His human variant. We demonstrate that the variant transcript is stably expressed and that agouti mice homozygote for the variant allele are viable with an altered coat color. This in vivo data confirms that the MFSD12 p.Tyr182His variant functions as a hypomorphic allele sufficient to alter mammalian pigmentation.

Lethal Yellow Mutation Causes Anxiety, Obsessive-compulsive Behavior and Affects the Brain Melanocortin System in Males and Females of Mice.

BACKGROUND: The brain melanocortin system regulates numerous physiological functions and kinds of behavior. The agouti protein inhibits melanocortin receptors in melanocytes. The lethal yellow (AY) mutation puts the Agouti gene under the control of the Raly gene promotor and causes the agouti protein expression in the brain. In the present article, we investigated the effects of the AY mutation on brain mRNA levels of Agouti, Raly, and melanocortin-related genes such as Agrp, Pomc, Mc3r, Mc4r, and their relationship to behavior. METHODS: The experiment was performed on 6-month-old males and females of AY/a and a/a (control) mice. Anxiety and obsessive-compulsive behavior were studied in elevated plus-maze and marble- burying tests. The mRNA levels were quantified by qPCR. RESULTS: AY mutation caused anxiety in males and obsessive-compulsive behavior in females. Positive correlation between Agouti and Raly genes mRNA levels were shown in the hypothalamus, hippocampus, and frontal cortex in AY/a mice. Reduced RNA concentrations of Mc3r and Mc4r genes were found respectively in the hypothalamus and frontal cortex in AY/a males. The Raly gene expression positively correlates with mRNA concentrations of the Mc3r gene in the hypothalamus and the Mc4r gene in the hypothalamus and frontal cortex. CONCLUSION: Possible association of obsessive-compulsive behavior with reduced Raly, Mc3r, or Mc4r gene expression is suggested.

The evolution of white-tailed jackrabbit camouflage in response to past and future seasonal climates.

The genetic basis of adaptive traits has rarely been used to predict future vulnerability of populations to climate change. We show that light versus dark seasonal pelage in white-tailed jackrabbits (Lepus townsendii) tracks snow cover and is primarily determined by genetic variation at endothelin receptor type B (EDNRB), corin serine peptidase (CORIN), and agouti signaling protein (ASIP). Winter color variation was associated with deeply divergent alleles at these genes, reflecting selection on both ancestral and introgressed variation. Forecasted reductions in snow cover are likely to induce widespread camouflage mismatch. However, simulated populations with variation for darker winter pelage are predicted to adapt rapidly, providing a trait-based genetic framework to facilitate evolutionary rescue. These discoveries demonstrate how the genetic basis of climate change adaptation can inform conservation.

A Cre knockin mouse reveals specific expression of Agouti gene in mesenchymal lineage cells in multiple organs and provides a unique tool for conditional gene targeting.

The mouse Agouti gene encodes a paracrine signaling factor which promotes melanocytes to produce yellow instead of black pigment. It has been reported that Agouti mRNA is confined to the dermal papilla after birth in various mammalian species. In this study, we created and characterized a knockin mouse strain in which Cre recombinase was expressed in-frame with endogenous Agouti coding sequence. The Agouti-Cre mice were bred with reporter mice (Rosa26-tdTomato or Rosa26-ZsGreen) to trace the lineage of Agouti-expressing cells during development. In skin, the reporter was detected in some dermal fibroblasts at the embryonic stage and in all dermal fibroblasts postnatally. It was also expressed in all mesenchymal lineage cells in other organs/tissues, including eyes, tongue, muscle, intestine, adipose, prostate and testis. Interestingly, the reporter expression was excluded from epithelial cells in the above organs/tissues. In brain, the reporter was observed in the outermost meningeal fibroblasts. Our work helps to illustrate the Agouti expression pattern during development and provides a valuable mouse strain for conditional gene targeting in mesenchymal lineage cells in multiple organs.

Aberrant expression of agouti signaling protein (ASIP) as a cause of monogenic severe childhood obesity.

Here we report a heterozygous tandem duplication at the ASIP (agouti signaling protein) gene locus causing ubiquitous, ectopic ASIP expression in a female patient with extreme childhood obesity. The mutation places ASIP under control of the ubiquitously active itchy E3 ubiquitin protein ligase promoter, driving the generation of ASIP in patient-derived native and induced pluripotent stem cells for all germ layers and hypothalamic-like neurons. The patient's phenotype of early-onset obesity, overgrowth, red hair and hyperinsulinemia is concordant with that of mutant mice ubiquitously expressing the homolog nonagouti. ASIP represses melanocyte-stimulating hormone-mediated activation as a melanocortin receptor antagonist, which might affect eating behavior, energy expenditure, adipocyte differentiation and pigmentation, as observed in the index patient. As the type of mutation escapes standard genetic screening algorithms, we rescreened the Leipzig Childhood Obesity cohort of 1,745 patients and identified four additional patients with the identical mutation, ectopic ASIP expression and a similar phenotype. Taken together, our data indicate that ubiquitous ectopic ASIP expression is likely a monogenic cause of human obesity.

Whole-Genome Sequencing Reveals the Genomic Characteristics and Selection Signatures of Hainan Black Goat.

Goats have become one of the most adaptive and important livestock species distributed in developing countries in recent years. The Hainan Black goat is a native goat breed of the Hainan region that is generally well-liked by the local population and is thus raised in large numbers. However, the genomic diversity and selective signals of the Hainan Black goat have not been clearly elucidated yet. Therefore, in this study, we performed whole-genome resequencing of 16 Hainan Black goats and compared the results with those of 71 goats of 6 other breeds from different geographic regions. Principal component analysis (PCA) and phylogenetic analysis identified seven lineages for all goats. Hainan Black goats showed the most similarity with Leizhou goats and the least similarity with Boer goats. Selective sweep analysis identified candidate genes associated with various functions, including immune resistance to disease (TNFAIP2 (TNF alpha induced protein 2) and EXOC3L4 (exocyst complex component 3 like 4)), melanin biosynthetic process (CDH15 (cadherin 15), ASIP (agouti signaling protein), and PARD3 (par-3 family cell polarity regulator)), and light sensitivity (CNGB3 (cyclic nucleotide gated channel subunit beta 3) and CNBD1 (cyclic nucleotide binding domain containing 1)), underlying strong selection signatures in Hainan Black goats. The melanin biosynthetic process, circadian entrainment, regulation of cyclic adenosine 3,5-monophosphate (cAMP)-mediated signaling, and the Rap-1 signaling pathway were significantly enriched in Hainan Black and Alashan Cashmere goats. This result may be important for understanding each trait. Selection signature analysis revealed candidate single nucleotide polymorphisms (SNPs) and genes correlated with the traits of Hainan Black goats. Collectively, our results provide valuable insights into the genetic basis of specific traits correlated with the Hainan island climate, artificial selection in certain local goat breeds, and the importance of protecting breed resources.

An enhancer of Agouti contributes to parallel evolution of cryptically colored beach mice.

Identifying the genetic basis of repeatedly evolved traits provides a way to reconstruct their evolutionary history and ultimately investigate the predictability of evolution. Here, we focus on the oldfield mouse (Peromyscus polionotus), which occurs in the southeastern United States, where it exhibits considerable color variation. Dorsal coats range from dark brown in mainland mice to near white in mice inhabiting sandy beaches; this light pelage has evolved independently on Florida's Gulf and Atlantic coasts as camouflage from predators. To facilitate genomic analyses, we first generated a chromosome-level genome assembly of Peromyscus polionotus subgriseus. Next, in a uniquely variable mainland population (Peromyscus polionotus albifrons), we scored 23 pigment traits and performed targeted resequencing in 168 mice. We find that pigment variation is strongly associated with an ∼2-kb region ∼5 kb upstream of the Agouti signaling protein coding region. Using a reporter-gene assay, we demonstrate that this regulatory region contains an enhancer that drives expression in the dermis of mouse embryos during the establishment of pigment prepatterns. Moreover, extended tracts of homozygosity in this Agouti region indicate that the light allele experienced recent and strong positive selection. Notably, this same light allele appears fixed in both Gulf and Atlantic coast beach mice, despite these populations being separated by >1,000 km. Together, our results suggest that this identified Agouti enhancer allele has been maintained in mainland populations as standing genetic variation and from there, has spread to and been selected in two independent beach mouse lineages, thereby facilitating their rapid and parallel evolution.

A genome-wide epistatic network underlies the molecular architecture of continuous color variation of body extremities.

Deciphering the molecular architecture of coat coloration for a better understanding of the biological mechanisms underlying pigmentation still remains a challenge. We took advantage of a rabbit French experimental population in which both a pattern and a gradient of coloration from white to brown segregated within the himalayan phenotype. The whole experimental design was genotyped using the high density Affymetrix® AxiomOrcun™ SNP Array and phenotyped into 6 different groups ordered from the lighter to the darker. Genome-wide association analyses pinpointed an oligogenic determinism, under recessive and additive inheritance, involving genes already known in melanogenesis (ASIP, KIT, MC1R, TYR), and likely processed pseudogenes linked to ribosomal function, RPS20 and RPS14. We also identified (i) gene-gene interactions through ASIP:MC1R affecting light cream/beige phenotypes while KIT:RPS responsible of dark chocolate/brown colors and (ii) a genome-wide epistatic network involving several others coloration genes such as POT1 or HPS5. Finally, we determined the recessive inheritance of the English spotting phenotype likely involving a copy number variation affecting at least the end of the coding sequence of the KIT gene. Our analyses of coloration as a continuous trait allowed us to go beyond much of the established knowledge through the detection of additional genes and gene-gene interactions that may contribute to the molecular architecture of the coloration phenotype.

Phylogeographic study using nuclear genome sequences of Asip to infer the origins of ventral fur color variation in the house mouse Mus musculus.

While the house mouse (Mus musculus), widely distributed in Eurasia, is known to have substantial coat color variation between and within local populations, in both primary and secondary distribution areas, including the Japanese archipelago, the evolutionary history of the color variation is poorly understood. To address the ventral fur color variation, we quantified the lightness of museum skin specimens, and found that the southern subspecies, M. m. castaneus (CAS), has high and low lightness in dry and rainy geographic regions, respectively. The northern subspecies, M. m. musculus (MUS), has low and high levels of lightness in the high and middle latitudes of northern Eurasia, respectively. We examined sequence variation of the agouti signaling protein gene (Asip), which is known to be responsible for the ventral fur color. We performed phylogenetic analyses with 196 haplotype sequences of Asip (~180 kb) generated by phasing the whole-genome data of 98 wild mice reported previously. Network and phylogenetic tree construction revealed clustering of haplotypes representing the two subspecies, MUS and CAS. A number of subclusters with geographic affinities appeared within the subspecies clusters, in which the essential results were consistent with those reconstructed with whole mitochondrial genome data, indicating that the phased haplotype genome sequences of the nuclear genome can be a useful tool for tracing the dispersal of geographical lineages. The results of phylogeographic analysis showed that CAS mice with darker ventral fur possessed similar Asip haplotypes across the geographic distribution, suggesting that these haplotypes are major causes of the historical introduction of Asip haplotypes for darker ventral fur in mice from northern India to the peripheral areas, including the Japanese archipelago. Similarly, MUS in East Asia, which has a white abdomen, formed an Asip haplogroup with that from northern Iran, also with a white abdomen.

The ASIP gene in the llama (Lama glama): Alternative transcripts, expression and relation with color phenotypes.

The Agouti gene (ASIP) is one of the most important genes for coat color determination in mammals. It has a complex structure with several promoters and alternative non-coding first exons that are transcribed into mRNAs with different 5'UTR. These mRNA isoforms regulate the temporal and spatial expression of the gene, producing diverse pigmentation patterns. Here, we studied ASIP transcriptional variants and their expression in the skin of llamas with different coat color phenotypes. We also described the ASIP locus, including promoter usage and the splicing events that originate each transcript variant. Using 5'RACE-PCR we isolated seven ASIP transcripts with alternative 5'UTR, where exons 1A, 1A', 1C, 1D, and a novel non-coding exon 1A" were identified. Additionally, new alternative spliced forms were found. The diversity of ASIP 5'UTRs is originated by a complex pattern of alternative promoter usage, multiple transcription start sites and splicing events that include exon skipping and alternative 3' splicing site selection. We found that ASIP was highly expressed in llamas with white and brown phenotypes while black animals presented very low expression. The main responsible for this difference was a fusion transcript between ASIP and NCOA6 genes, which was present in the skin of white and brown llamas but not in the black ones. The rest of ASIP transcripts presented very low expression in the skin, indicating that the main regulation point for ASIP gene expression is at the transcriptional level. Nevertheless, the characteristics of the 5'UTRs sequences suggest that alternative transcripts could be regulated differently at the protein synthesis level.

Cis-acting mutation affecting GJA5 transcription is underlying the Melanotic within-feather pigmentation pattern in chickens.

Melanotic (Ml) is a mutation in chickens that extends black (eumelanin) pigmentation in normally brown or red (pheomelanin) areas, thus affecting multiple within-feather patterns [J. W. Moore, J. R. Smyth Jr, J. Hered. 62, 215-219 (1971)]. In the present study, linkage mapping using a back-cross between Dark Cornish (Ml/Ml) and Partridge Plymouth Rock (ml+/ml+ ) chickens assigned Ml to an 820-kb region on chromosome 1. Identity-by-descent mapping, via whole-genome sequencing and diagnostic tests using a diverse set of chickens, refined the localization to the genomic region harboring GJA5 encoding gap-junction protein 5 (alias connexin 40) previously associated with pigmentation patterns in zebrafish. An insertion/deletion polymorphism located in the vicinity of the GJA5 promoter region was identified as the candidate causal mutation. Four different GJA5 transcripts were found to be expressed in feather follicles and at least two showed differential expression between genotypes. The results showed that Melanotic constitutes a cis-acting regulatory mutation affecting GJA5 expression. A recent study established the melanocortin-1 receptor (MC1R) locus and the interaction between the MC1R receptor and its antagonist agouti-signaling protein as the primary mechanism underlying variation in within-feather pigmentation patterns in chickens. The present study advances understanding the mechanisms underlying variation in plumage color in birds because it demonstrates that the activity of connexin 40/GJA5 can modulate the periodic pigmentation patterns within individual feathers.

Factors affecting gene expression associated with the skin color of black-bone chicken in Thailand.

The objective of this study was to investigate the effect of breed, sex, and age on the gene expression level of melanocortin 1 receptor (MC1R), DOPA chrome tautomerase (DCT), tyrosinase-related protein 1 (TYRP1), tyrosinase (TYR), and agouti signaling protein (ASIP) genes in Thai commercial chicken lines. All chicken have received Newscastle vaccination, and no antibiotics or any drugs were used in this study. Four chicken breeds including Black-Chinese, KU-Phuparn, Sri Mok, and Pradu Hang Dam were used in this study. These breeds can be classified by their skin color into 3 group including black (Black Chinese and KU-Phuparn), light black (Sri Mok), and yellowish white (Pradu Hang Dam). One hundred chickens per breed were used in this study. Breast skin tissue was randomly collected from 8 chickens (4 males, 4 females) per breed at 4, 8, 12, and 16 wk of age. The mRNA expression was analyzed using qRT-PCR and the gene expression level was calculated as 2-ΔΔCT. From the results, breed significantly (P < 0.01) affected the expression level for the 5 genes evaluated. Birds with the black skin color had greater TYRP1 and TYR gene expression when compared to chickens with light black and yellowish-white skin color, respectively. Whereas, chickens with yellowish-white skin color had greater ASIP gene expression when compared to chickens having the other skin colors. Sex significantly affected DCT, TYRP1, and TYR gene expression where the gene expression in males was greater when compared to females (P < 0.05). Age affected all gene expression levels (P < 0.01). At 4 wk of age, MC1R, DCT, TYRP1, and TYR gene expression was the highest and decreased as bird age increased (P < 0.05); however, ASIP gene expression was greatest at 8 wk of age. After 8 wk of age all gene expression for the genes evaluated in this study decreased as age increased. In addition, an interaction between breed and sex (P < 0.05) impacted DCT and ASIP gene expression. The results from this study showed that all genes evaluated can be used as candidate markers to further improve the blackness of the chicken's skin because the most desired skin color is black in the Thai black-bone chicken population.

Variants at the ASIP locus contribute to coat color darkening in Nellore cattle.

BACKGROUND: Nellore cattle (Bos indicus) are well-known for their adaptation to warm and humid environments. Hair length and coat color may impact heat tolerance. The Nellore breed has been strongly selected for white coat, but bulls generally exhibit darker hair ranging from light grey to black on the head, neck, hump, and knees. Given the potential contribution of coat color variation to the adaptation of cattle populations to tropical and sub-tropical environments, our aim was to map positional and functional candidate genetic variants associated with darkness of hair coat (DHC) in Nellore bulls. RESULTS: We performed a genome-wide association study (GWAS) for DHC using data from 432 Nellore bulls that were genotyped for more than 777 k single nucleotide polymorphism (SNP) markers. A single major association signal was detected in the vicinity of the agouti signaling protein gene (ASIP). The analysis of whole-genome sequence (WGS) data from 21 bulls revealed functional variants that are associated with DHC, including a structural rearrangement involving ASIP (ASIP-SV1). We further characterized this structural variant using Oxford Nanopore sequencing data from 13 Australian Brahman heifers, which share ancestry with Nellore cattle; we found that this variant originates from a 1155-bp deletion followed by an insertion of a transposable element of more than 150 bp that may impact the recruitment of ASIP non-coding exons. CONCLUSIONS: Our results indicate that the variant ASIP sequence causes darker coat pigmentation on specific parts of the body, most likely through a decreased expression of ASIP and consequently an increased production of eumelanin.

Functional conservation and divergence of color-pattern-related agouti family genes in teleost fishes.

While color patterns are highly diverse across the animal kingdom, certain patterns such as countershading and stripe patterns have evolved repeatedly. Across vertebrates, agouti-signaling genes have been associated with the evolution of both patterns. Here we study the functional conservation and divergence by investigating the expression patterns of the two color-pattern-related agouti-signaling genes, agouti-signaling protein 1 (asip1) and agouti-signaling protein 2b (asip2b, also known as agrp2) in Teleostei. We show that the dorsoventral expression profile of asip1 and the role of the "stripe repressor" asip2b are shared across multiple teleost lineages and uncover a previously unknown association between stripe-interstripe patterning and both asip1 and asip2b expression. In some species, including the zebrafish (Danio rerio), these two genes show complementary and overlapping expression patterns in line with functional redundancy. Our results thus suggest how conserved and novel functions of agouti-signaling genes might have shaped the evolution of color patterns across teleost fishes.

Neuropeptide Y-Agouti related peptide ratio (NAR) in patients with idiopathic primary hypothyroidism: Nudge and Risk.

OBJECTIVE: To illustrate the role of neuropeptide Y (NPY), agouti-related-peptide (AgRp) and their ratio in patients with primary hypothyroidism (PHT) regarding the effect of levothyroxine (LT4) replacement therapy. METHODS: The case-control study was conducted at the Department of Pharmacology, College of Medicine, Mustansiriya University, Baghdad, Iraq, from March to June, 2020 and involved 40 patients with primary hypothyroidism (PHT), including 20 newly diagnosed patients and 20 patients on levothyroxine (LT4) replacement therapy compared to 20 healthy controls. Anthropometric, lipid and pressure profiles were evaluated. Also T3, T4, TSH, NPY and AgRp serum level were estimated in different treated groups. SPSS version 20.00 was used for data analysis. RESULTS: Body mass index (BMI) was higher in the newly diagnosed patients without thyroxine therapy as compared to patients on thyroxine therapy (P=0.03). Blood pressure profile was higher in patients with PHT compared to the controls (P=0.0001). NPY serum level was lower in patients with PHT without thyroxine therapy (27.32±10.30ng/dL) as compared to patients with PHT on thyroxine therapy (61.10±22.78ng/dL), (P=0.04). AgRP serum level was lower in patients with PHT (9.81±4.86ng/dL) as compared to the patients on the thyroxine therapy (28.99±2.16ng/dL), (P=0.03). Besides, NPY-AgRP ratio (NAR) was higher in patients with PHT (2.78±0.14) as compared to patients with PHT on thyroxine therapy (2.10±0.19), (P=0.0001). CONCLUSION: Both of NPY and AgRP serum levels are reduced in the newly diagnosed patients with PHT and ameliorated following LT4 replacement therapy. Also NPY/AgRP ratio is linked with early PHT and regarded as a prognostic value for the outcomes of patients with PHT.

Investigating DNA methylation as a potential mediator between pigmentation genes, pigmentary traits and skin cancer.

Pigmentation characteristics are well-known risk factors for skin cancer. Polymorphisms in pigmentation genes have been associated with these traits and with the risk of malignancy. However, the functional relationship between genetic variation and disease is still unclear. This study aims to assess whether pigmentation SNPs are associated with pigmentary traits and skin cancer via DNA methylation (DNAm). Using a meta-GWAS of whole-blood DNAm from 36 European cohorts (N = 27,750; the Genetics of DNA Methylation Consortium, GoDMC), we found that 19 out of 27 SNPs in 10 pigmentation genes were associated with 391 DNAm sites across 30 genomic regions. We examined the effect of 25 selected DNAm sites on pigmentation traits, sun exposure phenotypes and skin cancer and on gene expression in whole blood. We uncovered an association of DNAm site cg07402062 with red hair in the Avon Longitudinal Study of Parents and Children (ALSPAC). We also found that the expression of ASIP and CDK10 was associated with hair colour, melanoma and basal cell carcinoma. Our results indicate that DNAm and expression of pigmentation genes may play a role as potential mediators of the relationship between genetic variants, pigmentation phenotypes and skin cancer and thus deserve further scrutiny.